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DETERMINATION  OF  VORTEX  PATHS  BY  SERIES  EXPANSION 
TECHNIQUE  WITH  APPLICATION  TO  CRUCIFORM  WINGS 
By  ALberta  Y.  Alksne 

SUMMARY 


A  series  method  of  determining  two-dimensional  vortex  paths  is  con¬ 
sidered  and  applied  to  the  computation  of  vortex  positions  behind  a  slender 
equal-span  cruciform  wing  at  any  angle  of  hank  as  a  function  of  the  dis¬ 
tance  behind  the  trailing  edge.  Calculated  paths  are  shown  for  four  hank 
angles.  For  a  hank  angle  of  45°  comparison  is  made  with  the  results  of  a 
closed  expression  given  in  NACA  TN  2605.  For  other  hank  angles  water-tank 
experiments  provide  qualitative  comparison.  Satisfactory  agreement  is 
found  for  a  sufficient  distance  downstream  to  include  most  practical 
missile -tail  positions. 

The  interference  forces  on  an  equal-span  interdigitated  cruciform  tail 
behind  a  slender  equal-span  cruciform  wing  are  calculated  for  four  angles 
of  hank  from  the  vortex  positions  found  by  use  of  the  series. 


INTRODUCTION 

It  is  now  well  established  that  the  vortex  wake  at  the  tail  of  a 
slender  configuration  similar  to  those  used  for  many  missiles  is  often 
entirely  rolled  up  and  that  the  downwash  field  at  the  tail  can  he  obtained 
by  use  of  a  single  discrete  vortex  as  an  approximation  to  the  vortex  wake 
trailing  from  each  wing  panel.  If  attention  is  confined  to  configurations 
for  which  the  ideas  of  conventional  slender -body  theory  can  be  used,  the 
problem  of  determining  the  steady-state  vortex  paths  becomes  an  exact 
analog  of  the  classical  problem  of  the  motion  of  a  system  of  parallel 
rectilinear  vortices . 

Sacks,  in  reference  1,  has  investigated  the  case  of  an  equal-span 
cruciform  wing  at  45°  angle  of  bank  where  the  symmetry  of  the  problem  per¬ 
mits  the  writing  of  a  closed  analytic  solution  for  the  vortex  paths.  The 
direct  extension  of  his  method  to  other  angles  of  bank  where  no  such 
symmetry  exists  does  not  appear  feasible. 

In  the  present  paper,  in  order  to  avoid  the  requirement  of  symmetry, 
a  series  has  been  developed  to  define  the  vortex  paths.  Paths  computed  by 
this  method  are  compared  with  the  analytic  results  of  Sacks  for  45°  angle 
of  bank,  and  with  the  results  of  water-tank  experiments  for  three  other 
bank  angles.  Calculations  are  made  of  the  forces  on  a  tail  due  to  vortices 
in  the  computed  positions . 
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SYMBOLS 

aspect  ratio,  — — 

S 

one  half  the  distance,  at  t  =  0,  between  the  two  vortices 
associated  with  a  component  wing1 

distance  of  a  vortex  from  the  plane  of  symmetry  of  a  component 
tail1 

maximum  chord 


lift  coefficient, 


(The  reference  area  used  in  this  report  is  the  area  of  one 
component  wing.) 

interference  lift  coefficient  (approximate) 

interference  normal-force  coefficient 

tail  surface  that  is  horizontal  at  0=5-  for  interdigitated 
tail  4 

perpendicular  distance  from  a  vortex  to  a  component  tail 

lift,  force  in  the  z  direction 

force  in  the  z  direction  on  a  cruciform  wing  (invariant  with 
bank  angle) 

projection  of  the  interference  normal  force  on  the  x,z  plane 
(approximately  the  interference  lift) 

interference  normal  force,  that  is,  normal  force  on  a  component 
tail  due  to  the  presence  of  vortices 

positive  integer 

radius  of  a  cylindrical  boundary;  specifically,  the  radius  of  the 
water  tank 

•JJy7  -  yJ2  +  (z,  -  z.)2 


(The  reference  area  used  in  this  report  is  the  area  of  one 
component  wing . ) 


xBy  "component  wing"  (or  "component  tail")  is  meant  a  wing  (or  tail) 
consisting  of  two  "panels"  lying  in  the  same  plane. 
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ANALYSIS 
Axis  System 

The  coordinate  system  used  in  this  report  is  a  wind-axis  system  as 
shown  in  figure  1,  where  the  origin  lies  at  the  trailing  edge  of  the  wing 
center  line.  The  angle  of  attack  is  required  hy  the  limitation  of  the 
theory  to  he  small  and  the  trailing  edges  of  the  wing  panels  are  assumed 
to  lie  in  the  x  =  0  plane. 


Series  Solution  for  Vortex  Motions 


Analytical  solutions  for  the  motion  of  a  system  of  parallel  recti^ 
linear  vortices  are  given  hy  Grohli  in  reference  2  for  the  case  of  three 
vortices  with  certain  restrictions  on  the  starting  positions  and  strengths 
for  four  vortices  with  a  plane  of  symmetry,  and  for  2n  vortices  with  n 
planes  of  symmetry.  The  solution  for  four  vortices  as  given  hy  Grohli 
contains  an  errors  hut  is  given  correctly  hy  Sacks  in  reference  1  and  is 
there  applied  to  the  case  of  the  vortices  behind  a  slender  equal-span 
cruciform  wing  at  ^5°  angle  of  hank,  that  is,  to  four  vortices  of  equal 
strength  starting  in  the  form  of  a  square.  The  solution  in  this  ‘case 
depends  on  the  existence  of  a  plane  of  symmetry  and  cannot  readily  he 
extended  to  cases  of  arbitrary  vortex  strength  where  the  symmetry  is 
lacking. 

The  present  analysis  undertakes  to  define  the  positions  of  a  number 
of  vortices  of  given  strengths  and  initial  positions  in  terms  of  a  Taylor* 
series  in  powers  of  the  time,  t,  thereby  eliminating  the  dependence  on 
symmetry.  Expansion  around  t  =  0  results  in  the  following  expression 
for  the  position  of  the  i^h  vortex: 


The  coefficients  of  this  series  can  he  determined  hy  using  the  Biot- 

Savart  law  for  two-dimensional  vortices  parallel  to  the  x  axis.  For  a 

system  of  free  vortices,  if  v-j_  and  Wj[  are  the  y  and  z  components  of 

the  velocity  of  the  i*th  vortex  due  to  a  vortex  of  strength  Tj  situated 

at  y^,  z^,  the  required  vortex  laws  are: 
d  7  d 

2Page  1^7  of  reference  2,  equations  23  and  2^. 
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vi(t)  =  T7  yi(t)  -  > 


r n  z^( t)  -  zi(t) 


w  m  d  „  m  V  r.i  yi(t)  "  yJ(t) 

wi'^'  ~  TT  zivt'  ~  )  o«  ~p 


where 


■2  =  [y±(t) 


yj(t)]2  +  [z±(t)  -  Zj(t)] 


Now  if  the  positions  of  all  the  vortices  are  known  at  t  =  0,  it  is 
possible  to  write  the  coefficient  of  the  first  power  of  t  in  equation  (l) 
for  all  the  vortices  concerned  by  simply  substituting  the  initial  positions 


clz 

into  equation  (2)  to  get  (  — )  and  ( — i )  .  Thus  the  first  two  terms 

\dt/t=0  \  dt  / 1=0 

of  the  series  are  known  for  all  the  vortices,  that  is,  y^  and  Zj_  can  now 
be  written  as  linear  functions  of  t.  Substituting  these  first  two  terms 
into  equation  (2)  and  differentiating  with  respect  to  t  and  then  setting 
t  equal  to  zero  gives  the  coefficient  of  the  second  power  of  t.  Now 
three  terms  of  the  series  are  available  for  substitution  into  equation  (2), 
etc.  Note  that  at  each  step  the  unknown  terms  of  the  series- are  of  no 
significance  in  the  process  since  they  still  contain  t  as  a  factor  after 
the  differentiation  and  therefore  disappear  when  t  is  set  to  zero. 

The  following  formula  for  differentiation  of  a  product  of  two  func¬ 
tions  is  convenient  for  use  in  obtaining  higher  order  terms: 


Thus  the  first  two  terms 


™  g  *1  f  g  ♦ 

bn’x  21  dt2  dtn"2 


11;  a  a.  -p  1 

- g  -  f  +  •  . 

(n  -  k)  Ik'.  dtn“k  dtk 


.  g  —  f  (3) 
dtn 


n  >  1 


0  <  k  <  n 
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Now  it  can  be  seen  that  if  the  positions  can  actually  he  described 
by  such  a  series,  the  only  restriction  on  problems  to  be  solved  is  that 
the  series  should  converge  rapidly  enough  to  be  practical  for  the  desired 
values  of  t  and  that  the  work  of  evaluating  the  coefficients  should  not 
be  prohibitive. 

As  a  test  of  the  method  the  coefficients  have  been  determined  out  to 
the  fourth  power  of  t  for  the  case  which  corresponds  to  the  equal-span 
cruciform  wing  at  any  angle  of  bank,  that  is  for  four  vortices  initially 
placed  at  the  corners  of  a  square  with  diagonally  opposite  vortices  of 
equal  strength  but  opposite  sense.  Furthermore,  since  there  was  a  closed 
analytic  solution  available  for  this  configuration  at  an  angle  of  bank  of 
45  ,  five  additional  coefficients  were  found  for  that  case  with  a  view  to 
increased  understanding  of  the  behavior  of  the  series. 

Initial  Positions  and  Strengths  of  Vortices 

In  accordance  vith  the  work  of 
Spreiter  and  Sacks  (ref.  3)  all  of  the 
vorticity  behind  a  wing  has  been  assumed 
to  be  concentrated  in  vortex  lines 
springing  from  the  centroid-of -vorticity 
positions  at  the  trailing  edge  and  sub¬ 
ject  thereafter  to  the  two-dimensional 
vortex  laws.  Since  the  circulation,  r, 
y  is  equal  to  the  jump  in  potential,  Acp, 
and  since  slender -body  theory  leads  to 
an  elliptic  spanwise  distribution  of 
Acp,  the  centroid  of  vorticity  at  the 
trailing  edge  of  each  wing  panel  lies 
at  a  point  n/k  of  the  distance  from 
the  center  line  to  the  wing  tip.  Thus, 
each  vortex  pair  has  a  span  of  2a,  where 
a  =  (jt/U)sw  as  shown  in  sketch  (a). 


Sketch  (a) 
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The  lift  of  the  equal-span  cruciform  wing,  1^,  that  is,  the  force 
in  the  direction  of  the  positive  z  axis  does  not  vary  with  angle  of  hank 
hut  remains  throughout: 


Lw  =  PjyV^a) 


(4) 


where  rw  is  the  circulation  around  the  horizontal  wing  at  angle  of  hank 
0  =  0.  At  other  angles  of  hank  the  vortex  strengths  are  related  to  Tw 
as  follows: 


r4  =  "  r2  =  rwcos  0 

ra  =  -  r3  =  rvsin  0 


(5) 


where  the  vortices  are  numbered  as  in  sketch  (a). 


Solutions 


The  use  of  equations  (l),  (2),  (3),  and  (5),  together  with  the  fact 
that  the  initial  positions  of  the  vortices  are  known  in  terms  of  the  wing 
semispan,  sw,  and  the  hank  angle,  0,  leads  to  a  series  for  the  vortex 
positions  at  any  time  t.  Coefficients  have  been  found  out  to  the  fourth 
power  of  t,  and  with  the  substitution 


rw* 

kna2 


(6) 


the  series  can  he  written  as  follows: 
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liy 

— -  =  sin  0  -  (sin  0  cos  0)T  -  £  (sin  0  cos30)T3  + 

jrsv  3 

—  (sin  0  cos20)(sin20  -  cos20)T4.  .  . 


=  cos  0  -  (l  +  cos20)T  -  (cos  0)T2  -  i  cos20(l  +  2  cos20)T3 
i  cos  0(3  -  4  sin20  cos20)T4.  .  . 


izk 


4z, 


-  cos  0  +  (sin  0  cos  0)T  +  ^  (sin30  cos  0)T3  + 
~  (sin20  cos  0)(sin20  -  cos20)T4.  .  . 


— ^  =  sin  0  -  (l  +  sin20)T  -  (sin  0)T2 
*s v 

~  sin  0(3-4  sin20  cos20)T4.  .  . 


—  sin20(l  +  2  sin20)T3 


4y3 

rtsW 


4z 


-  sin  0  -  (sin  0  cos  0)T  -  2  (sin  0  cos30)T3 
~  (sin  0  cos20)(sin20  -  cos20)T4.  .  . 


>  (7) 


=  -  cos  0  -  (l  +  cos20)T  +  (cos  0)T2  -  i  cos20(l  +  2  cos20)T3  + 
i  cos  0(3-4  sin20  cos20)T4.  •  . 


=  cos  0  +  (sin  0  cos  0)T  +  —  (sin30  cos  0)T3  - 
itsw  3 

^  (sin20  cos  0)(sin20  -  cos20)T4.  .  . 


4z4 

JtSw 


=  -  sin  0  -  (l  +  sin20)T  +  (sin  0)T2  -  i  sin20(l  +  2  sin20)T3  + 
i  sin  0(3-4  sin20  cos20)T4.  .  • 
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For  the  special  case  treated  by  Sacks,  that  is,  the  equal-span 
cruciform  wing  at  an  angle  of  bank  0  =  jt/4,  there  is  a  plane  of  symmetry 
and  the  series  is  simplified  so  that  the  labor  of  evaluating  the  coeffi¬ 
cients  is  a  great  deal  less.  For  this  case  coefficients  were  found  out 
to  the  ninth  power  of  t.  Since  cos  0  =  sin  0  =  l/s/2,  and 


r t  =  -  r,  «  r,  = 


F 2  =  the  series  can  be  written,  letting  r  =  t//2: 

V2 


kj2Y1_  Q  3 

*  sw  3 


A  T5  +  iA  tt  +  13318  9 
15  315  2835 


i^^=l-3r 

It 


2t2  -  i  T3  -  8  t4  _  28  Ts 

3  3  15 

3062  „  12172  „ 

-  t8  -  - - -  T9  •  •  . 

315  2835 


226  s 

45 


itil  Zi  =  1  +  T  *  2  T3  t  i.  T=  .  13i  t7  .  13215  tS.  . 

It  Sv  3  15  315  2835 


4  -v/2  Z4 

K  Sv 


-1  -  3t  +  2r2  -  -  T3  +  -  T4  -  —  T5  + 
3  3  15 

t6  .  932  t7  +  3062  t8  _  12172  tS> 
45  315  315  2835 


Relation  Between  T  and  x 


As  pointed  out  previously,  the  vortex  laws  used  here  apply  to 
straight-line  vortices,  parallel  to  the  x  axis,  extending  to  infinity 
in  both  directions,  and  changing  their  position  with  time.  Their  use  in 
the  situation  to  which  they  are  to  be  applied  depends  on  the  three- 
dimensional  steady-state  vortex  picture  showing  relatively  gradual  vari¬ 
ations  in  the  x  direction.  Within  the  limits  of  slender -body  theory 
the  correspondence  is  exact,  and  the  results  obtained  in  the  previous 
section  can  be  used  to  compute  the  three-dimensional  vortex  paths  behind 
a  slender  equal-span  cruciform  wing  by  means  of  the  relation 


x  =  ^t 


(9) 
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Now  since 

T  =  rvt 

krta2 


(6) 


and 


Lw  =  PooUoorv( 2a) 


W 


and  in  slender  theory 


=  (Tt/2)Aa  (10) 

the  positions  of  the  vortices  at  any  downstream  station,  x/sv,  can  he 
found  from  equation  (7),  or  for  45°  hank  from  equation  (8),  hy  use  of 
the  relation 


x  _  *4(sw)2  T  _  ^(sy)2  T 
Sw  ,  2AySw^w 


where  5^  is  the  "attitude  angle,11  that  is,  the  angle  between  the  center 
line  of  the  cruciform  wing  and  the  free-stream  direction. 


If 

written: 


as  for  triangular  wings,  then  equation  (ll)  can  he 


—  =  T  =  -2^-  T  (12) 

sw  8^w 

Equation  (12)  provides  the  relation  used  in  the  present  report. 


EXPERIMENT 


In  order  to  provide  a  qualitative  means  of  judging  the  results  of 
the  computations  for  angles  of  hank  for  which  no  closed  analytic  solution 
is  available,  experiments  were  run  with  small  models  in  a  water  tank. 
Water  miscible  paint  spread  on  the  trailing  edge  before  each  run  remained 
floating  on  the  surface  of  the  water  behind  the  model  and  made  the  vor¬ 
tices  visible.  Eor  various  reasons,  it  was  considered  inadvisable  to 
attempt  quantitative  comparison.  For  one  thing,  there  is  no  general 
agreement  as  to  the  point  in  a  vortex  swirl  which  is  to  be  considered 
the  center  of  the  core,  and  the  centroid  of  vorticity,  which  is  the 
quantity  calculated  in  this  report,  is  even  more  difficult  to  define. 

For  another  thing,  the  best  pictures  were  obtained  at  angles  of  attack 
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which  were  too  high  to  he  entirely  compatible  with  the  assumptions  of 
the  theory.  However,  the  water-tank  experiments  were  expected  to  demon¬ 
strate  the  trends  in  the  variation  of  the  vortex  patterns  with  hank  angle 


Water  Tank  and  Models 


The  water  tank  used  in  the  present  experiments  was  the  same  as  that 
described  by  Sacks  (ref.  l)  and  is  shown  in  figure  2.  Three  different 
models  were  used,  all  equal-span  cruciform  wings  constructed  of  sheet 
metal  0.050  inch  thick.  One  model  had  an  8-inch  span  and  an  aspect  ratio 
of  2.  The  others  were  smaller,  having  only  a  4-inch  span.  One  of  these 
had  an  aspect  ratio  of  1,  the  other,  2.  Various  angles  of  attack  were 

tried.  The  most  successful  runs  were  made  at  a  =  jg  =  15° •  As  in  ref¬ 
erence  1,  motion  pictures  provided  a  record  of  the  distance  traveled  by 
the  wing  as  well  as  of  the  changing  vortex  patterns. 

The  water  tank  was  not  deep  enough  for  the  8-inch-span  model  to 
continue  running  much  beyond  two  span  lengths  below  the  surface .  However 
the  camera  was  kept  running  after  the  model  stopped  and  the  time  measured 
in  frames  was  used  to  determine  an  equivalent  distance . 


Accuracy  and  Repeatability  of  Experimental  Data 


It  was  found  that  runs  made  with  the  two  small  models  showed  exces¬ 
sive  influence  of  currents  set  up  in  the  tank  by  the  supporting  mechanism 
and  by  various  outside  distrubances .  The  vortex  paths  behind  the  model 
with  8-inch  span  showed  little  effect  of  stray  disturbances  but,  since 
the  water -tank  diameter  was  only  22  inches,  there  was  a  large  systematic 
error  due  to  wall  interference. 

A  comparison  was  made  of  the  vortex  patterns  at  corresponding  dis¬ 
tances  behind  the  three  different  models  for  0  =  jt/4.  The  choosing  of 
the  particular  runs  in  which  the  symmetrical  vortex  pattern  typical  of 
this  symmetrical  configuration  was  maintained  eliminated  most  of  the 
irrelevant  disturbances  which  made  data  from  the  small  models  generally 
unsatisfactory.  Measured  in  terms  of  the  half  span  of  the  model,  the  y 
coordinates  of  the  vortex  cores  at  corresponding  distances  behind  the 
three  models  did  not  differ  by  more  than  10  percent.  However,  the  z 
coordinates  behind  the  large  model  differed  from  those  behind  the  small 
models  by  about  25  percent. 

A  wall -interference  correction  consisting  of  a  constant  upwash, 
computed  on  the  assumption  of  four  discrete  vortices  (see  Appendix  A) 
was  sufficient  to  bring  the  results  for  the  large  model  into  very  good 
agreement  with  those  for  the  small  models.  Since  this  was  the  case  for 
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the  hank  angle  of  k^°  where  the  general  nature  of  the  vortex  pattern  was 
known,  and  since  the  upwash  at  t  =  0  was  found,  under  the  same  assump¬ 
tion,  not  to  vary  with  "bank  angle,  it  was  assumed  that  the  data  obtained 
from  the  large  model  for  other  angles  would  also  be  satisfactory  when  the 
same  correction  was  applied.  The  water -tank  pictures  shown  in  this  report 
are  those  taken  with  the  large-span  model  and  the  necessary  corrections 
are  indicated  by  additional  reference  points  marked  at  the  sides  of  the 
prints . 

No  allowance  was  made  for  the  effects  of  stopping  the  model  before 
the  runs  were  complete.  The  influence,  if  any,  should  have  appeared  as 
an  additional  downwash  at  the  surface  near  the  end  of  a  run,  but  none  was 
noted  in  comparing  runs  made  with  the  large  and  small  models. 


RESULTS  AND  DISCUSSION 


Computations  of  vortex  paths  behind  a  cruciform  wing  have  been  made 
using  equations  (7)  and  (8)  for  four  angles  of  bank,  0  =  ic/l 6,  jt/8,  3tc/1  6, 
jt/4.  Figure  3.  shows  the  paths  with  y/sv  plotted  against  z/sw  for 
various  values  of  x/sw  in  a  coordinate  system  in  which  the  x  axis 
lies  in  the  stream  direction  and  the  bank  angle  is  measured  from  the  z 
axis  and  is  positive  when  the  starboard  wing  is  rotated  down.  The  points 
shown  are  for  T  =  0,  0.276,  O.352,  0.449,  O.517,  O.582,  O.65O,  and  O.766, 
which  for  an  aspect-ratio-2  wing  with  a  lift  coefficient  of  0.82  (<x=jt/l2) 
corresponds  to  x/sw  =  0,  4.1,  5*2,  6.6,  7*6,  8.6,  9*6,  and  11.4;  that  is, 
this  figure  may  be  considered  either  as  a  time  history  or  as  a  representa¬ 
tion  of  three-dimensional  vortex  paths. 

Computations  have  been  carried  out  to  a  value  of  T  corresponding 
to  the  "leapfrog"  position  of  reference  1,  that  is,  to  the  value  at  which 
the  two  upper  vortices  pass  between  the  two  lower  vortices  for  a  bank 
angle  of  45°. 

Figure  4  shows  water-tank  pictures  taken  with  the  8-inch-span  model 
at  an  attitude  angle  of  tc/12  radians  and  at  bank  angles  of  jt/l6,  tc/8, 
and  3^/l6 .  Choice  of  the  appropriate  frames  from  the  motion  picture  film 
made  it  possible  to  present  pictures  corresponding  very  closely  to  most  of 
the  values  of  T  used  in  the  computations.  No  comparison  is  shown  for 
small  values  of  T  where  the  vortices  were  in  the  process  of  rolling  up 
and  the  visible  vortex  cores  were  not  only  poorly  defined  but  did  not  yet 
correspond  to  the  centroids  of  vorticity.  As  can  be  seen  from  the  first 
picture  of  each  series,  the  solid  white  markers  indicate  the  point  at 
which  the  trailing-edge  center  line  entered  the  water.  The  open  white 
markers  indicate  the  corrected  position  of  this  reference  point,  shifted 
upward  to  account  for  the  upwash  due  to  the  presence  of  the  tank  wall. 

(See  Appendix  A.)  It  can  be  seen  that  the  variation  with  angle  of  bank 
found  by  the  calculations  (fig.  3)  is  similar  to  that  shown  in  figure  4. 
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In  figure  3(d)  the  positions  as  calculated  by  the  formulas  of  ref¬ 
erence  1  are  shown  for  comparison  with  the  series  results  at  a  bank  angle 
of  jr/4.  For  this  bank  angle  the  ninth-order  terms  of  the  series  were 
available  from  equation  (8).  Points  are  also  shown  computed  with  terms 
out  to  the  fourth  order  as  for  the  other  angles  of  bank.  Even  at  the 
"leapfrog11  position  (last  point  computed)  the  agreement  is  good  if  ninth- 
order  terms  are  used.  If  only  terms  out  to  the  fourth  order  are  used, 
the  largest  error  appears  in  the  z  coordinate  of  the  fast  moving  vor¬ 
tices  1  and  2,  but  at  a  point  which,  for  5  =  12,  would  correspond  to  a 
distance  downstream  of  four  times  the  wing  span,  the  error  is  still  less 
than  10  percent  of  the  total  change  of  position  in  the  z  direction,  or 
about  5  percent  of  the  wing  span. 

From  this  comparison  with  the  work  of  Sacks,  together  with  the  fact 
that  the  water-tank  pictures  for  other  angles  of  bank  also'  show  qualita¬ 
tive  agreement  with  the  computed  vortex  positions,  it  appears  that  the 
series  computations  give  satisfactory  results  for  a  distance  of  several 
wing  spans  behind  the  trailing  edge. 


As  a  further  check  on  the  dependability  of  the  series  method  the 
results  for  0  =  ir/4  using  one  term  of  the  series,  two  terms,  three  terms 
etc.,  are  presented  in  figures  5  and  6,  again  in  comparison  with  the 
results  calculated  from  the  formulas 


of  Sacks.  In  these  figures  the 
values  of  y/sw  and  z/sw  are  plotted 
separately  against  the  parameter  T, 
which  is  related  to  the  downstream 
distance  as  in  equation  (12),  so 
that 

xCT  l6 


SwAwrt4 


It  can  be  seen  from  these  figures 
that  the  series  appears  to  converge 
quite  rapidly  for  small  values  of 
T,  and  to  converge,  although  more 
slowly,  even  for  the  highest  value 
of  T  used.  Sketch  (b)  shows  the 
sum  of  the  first  n  terms  plotted 
against  n,  out  to  n  =  9,  for 
T  =  0.517  for  vortices  1  and  4. 

(The  first  term,  n  =  0,  is  not  shown 
as  it  is  simply  the  initial  posi¬ 
tion.)  It  can  be  seen  that  in  both 
cases  the  series  for  y/sw  con¬ 
verges  very  quickly.  Note  also  that 
only  the  odd  powers  of  T  appear  in 


n 


4  =  tt/4|  T  =  0.517 


Sketch  (b) 
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the  series  for  y/sw.  The  series  for 
Z4/SV  alternates  and  that  for  zi/sw 
does  not;  however ,  both  approach  the 
correct  value  very  rapidly  at  T  =0.517* 

Figures  7  and  8  show  y/sw  and 
z/sw  plotted  against  T  for  each  of 
the  four  vortices  for  a  bank  angle  of 
tc/8.  Only  fourth -order  terms  are 
available  for  this  case,  but  it  can  be 
seen  that  the  behavior  of  the  series 
is  very  similar  to  that  observed  for 
0  =  rt/4.  Sketch  (c)  shows  the  sum  of 
the  first  n  terms  for  0  =  tc/8  plot¬ 
ted  against  n  out  to  n  =  4  for  each 
vortex  for  T  =  O.517,  further  sub¬ 
stantiating  the  statement  concerning 
similar  behavior,  although  the  con¬ 
vergence  is  slower. 


Sketch  (c) 

It  appears  from  these  considerations  that  for  small  values  of  T 
only  two  or  three  terms  of  the  series  are  required  and  that  four  terms 
are  enough  to  give  good  results  out  to  about  T  =  0.5.  Beyond  that  point 
the  results  become  somewhat  doubtful  unless  more  terms  are  used,  but  this 
provides  a  fairly  wide  range  of  useful  values,  as  can  be  seen  from  the 
fact  that  for  a  lift  coefficient  of  0.5  and  an  aspect  ratio  of  2,  T  =  0.5 
corresponds  to  about  six  span  lengths  behind  the  trailing  edge  of  the  wing. 


CALCULATION  OF  LIFT  ON  A  CRUCIFORM  TAIL 


The  lift  on  the  tail  of  a  slender  wing-tail  combination  due  to  the 
vortices  from  the  wing  can  be  computed  by  reverse  flow  techniques  as 
discussed  in  reference  3,  on  the  assumption  that  the  tail  does  not  influ¬ 
ence  the  positions  of  the  vortices.  The  equation 


%  = 


rP00U00st 
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given  in  reference  4  yields  the  nor¬ 
mal  force  on  a  component  tail,  or 
tail  plane  (see  sketch  (d)),  due  to 
a  single  vortex.  Note  that  a  factor 
has  teen  placed  In  front  of  the 
outer  radical  to  take  account  of  the 
case  of  t  <  0.  The  effects  of  all 
the  vortices  must  he  summed  for  each 
component  tail  and  the  components  of 
force  in  the  z  direction  added  to 
give  the  lift. 

The  configuration  chosen  for 
the  present  calculations  was  a  cru¬ 
ciform  tail  interdigitated  behind  a 
cruciform  wing  as  in  sketch  (e), 
where  the  tail  components  are  desig¬ 
nated  V  and  H  as  shown.  Note  that 
the  tail  center  line  is  an  extension 
of  the  wing  center  line,  whereas  the 
x  axis  lies  in  the  stream  direction. 
Since  no  account  was  taken  of  the 
effect  of  the  tail  on  the  vortex 
paths,  the  vortex  positions  used  in 
the  force  calculations  were  those 
already  computed  in  the  absence  of 
a  tail  for  the  station  corresponding 
to  the  tail  trailing  edge. 


z 


The  normal -force  coefficients  on  each  tail  plane  have  been  computed 
from  equation  (13)  for  four  angles  of  bank  for  three  ratios  of  tail  span 
to  wing  span,  and  for  T  =  0.247  and  T  =  0.4ll.  It  should  be  noticed 
that  varying  T  corresponds  to  varying  either  the  tail  length,  the  wing 
lift  coefficient,  or  the  aspect  ratio  (see  eq.  (12)).  For  Cl  =  0.5  and 
A  =  2  these  values  of  T  correspond  to  x/sw  =  6  and  x/sw  =  10.  The 
results  are  shown  in  figure  9* 


In  figure  10  is  shown  the  lift  coefficient  of  the  tail  due  to  the 
presence  of  the  vortices,  that  is,  the  interference  lift  coefficient, 
CLl,  measured  in  the  positive  z  direction  (the  same  direction  as  the 

lift  on  the  wing).  It  is  interesting  to  note  that  for  T  =  0.2^7  and 
st/sw  =  1,  and  also  for  T  =  0.4ll  and  st/sw  =  1.2,  there  is  very  little 
variation  of  interference  lift  coefficient  with  angle  of  hank. 

A  comparison  with  analytical  results  obtained  by  the  method  of  ref¬ 
erence  5  for  k-5°  angle  of  bank  is  shown  in  figure  11,  in  which  the  inter 
ference  lift  is  plotted  against  the  tail-span  to  wing-span  ratio.  The 
agreement  is  very  good,  as  of  course  should  be  expected  since  the  vortex 
positions  agree  so  well. 


1 6 
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CONCLUDING  REMARKS 


Computations  of  vortex  paths  and  of  forces  on  a  tail  behind  a  slender 
equal-span  cruciform  wing  have  been  made  using  a  series  to  determine  the 
vortex  positions  as  a  function  of  the  distance  downstream.  The  results 
show  that,  for  a  bank  angle  of  45° ,  only  a  few  terms  of  the  series  are 
needed  to  provide  satisfactory  agreement  with  the  known  analytic  solution 
at  downstream  distances  encompassing  most  practical  missile  tail  posi¬ 
tions.  Comparison  with  water -tank  pictures  of  the  vortex  patterns  for 
other  angles  of  bank,  and  consideration  of  the  relative  size  of  succes¬ 
sive  terms  of  the  series,  indicate  that  the  same  is  true  for  the  general 
case  where  no  analytic  solution  is  known. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Feb.  8,  195 6 
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APPENDIX  A 


WALL-INTERFERENCE  CORRECTIONS 


As  in  the  case  of  tunnel-wall  corrections ,  the  effect  of  surrounding 
a  group  of  vortices  with  a  solid  cylindrical  boundary  of  radius  R  can 
be  calculated  by  the  method  of  images  when  the  vortex  positions  and 
strengths  are  known. 

For  each  vortex  within  the  boundary,  the  position  is  given  as  £ . 
where  ^  =  y  +  iz  and  i  =  «/- 1 .  The  position  of  the  image  outside  the^ 
boundary  is  then  known  and  is  R2/f j  where  f  =  y  -  iz.  Then  at  a 
point  (;  the  complex  velocity  due  to  the  image  vortex  is 


(Al) 


If  the  field  of  interest  is  confined  to  a  small  area  in  the  center  of 
the  cylinder  and  if  the  vortices  also  remain  in  this  area,  £  may  be 
neglected  as  very  small  compared  to  R2/£j  and  the  complex  velocity 
due  to  a  number  of  image  vortices  may  be  written 


u  -  f  l(-rjL>  (A2) 

L>  2tiR2 

j=i 

Then 


n  r  _  n 

v  ■  X  S&  IP1?J  ■  t 

rjyj 

2«R2 

j=i 

n  (  n  .  1  n 

(-rjzj) 

2rtR2 

J=1  j=l 

(A3) 


Since  the  present  report  is  concerned  with  an  equal-span  cruciform 
wing,  there  are  assumed  to  be  only  four  vortices  within  the  boundary  and 
the  relation  between  their  strengths  is 
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ri  =  -r3  =  rwsin  $ 
r4  =  -r2  =  rwcos  0 


1 

(5) 


Then  at  the  center  of  the  cylinder 


v  =  £(yx  -  y3)sin  0  +  (y4  -  y2)cos 


V  = 


w 


2*R2 


(zi  -  z3)sin  0  +  (z4  -  z2)cos 


-] 


}  (A4) 


where  yx,  y2,  y3,  y4  and  zi,  z2,  z3,  and  z4  depend  on  t.  At  t  =  0 
the  positions  of  the  vortices  are  known  in  terms  of  0  and  the  semispan, 
sv,  and  the  expressions  for  w  and  v  simplify  to 


Fwsw 

UR2 


(A5) 


for  any  angle  of  hank. 

For  the  special  case  of  ^5°  hank  angle,  sin  0  =  cos  0  =  — and 
symmetry  provides  relations  between  the  vortex  positions  so  that 


w  =  —Ti=  (yi  +  y4) 

rtR2v2 


v  =  0 


(A6) 


t 


at  any  time  t.  Furthermore  it  is  known  (see  refs.  1  and  2)  that  for 
this  case  (y1  +  y4)  is  constant  with  time  so  at  the  center  of  the  cylinder 

w  =  constant;  v  =  0  (A7) 

as  long  as  symmetry  with  respect  to  the  z  axis  is  maintained. 
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Since,  near  the  center  of  the  cylinder,  the  upwash  for 
hank  is  the  same  at  t  =  0  and  the  upwash  for  0  =  Jt /k  does 
with  time,  it  has  been  assumed  that  one  correction,  namely. 


fwsw 

4r2 


could  he  used  throughout.  This  resulted  in  a  correction  of 
the  z  position  at  any  time,  t,  where  R  is  the  radius  of 
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the  water  tank. 
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Figure  2.-  Water  tank  with  cruciform  model 
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-1.0  -.5  0  .5  1.0 

y/sw 

(a)  0  =  jr/l6 

Figure  3.-  Vortex  paths  behind  a  cruciform  wing  at  four  different  angles 

of  hank . 
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(t>)  0  =  ir/8 
Figure  3*“  Continued. 
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(c)  0  =  3«/i6 
Figure  3.-  Continued. 
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(d)  0  =  n/h 
Figure  3,-  Concluded, 
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x/s^-  ~  ^ 


x/s w  =  4.1 
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X/Sy  =  7.6 


x/^  =  5.2  x/sw  =  8.6 

(a)  <f  =  )t/l6 

Figure  4.-  Photographs  of  wake  at  various  stations  behind  a  cruciform 
wing  of  aspect  ratio  2,  for  three  different  bank  angles;  a  =  tt/12. 


Figure  4.-  Continued, 
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x/sw  =9*6 

(t>)  0  =  it/8  -  Concluded. 


Figure  4.-  Continued. 


Figure  5.-  Variation  with  the  parameter  T  of  the  lateral  position  y/s- 
of  vortex  1  and  vortex  k  behind  a  cruciform  wing,  calculated  with 
one,  two,  three,  four,  etc.,  terms  of  the  series;  f  =  7t/4. 
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(a)  Vortex  number  I 
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(c)  Vortex  number  3 


Figure  8.-  Concluded, 
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interdigitated  cruciform  tail  at  two  values  of  the  downstream  distance  parameter  T  behind 


4o 
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Figure  11.-  Variation  with  the  ratio  of  tail  span  to  wing  span  of  the 
interference  lift  coefficient  on  an  interdigitated  cruciform  tail 
behind  a  cruciform  wing  at  angle  of  hank,  0  =  n/k;  CL  =  0.5;  x/sw  =  6. 
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